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We derive expressions for the intensity of the Brillouin polarized spectrum of a molecular liquid 
formed of axially symmetric molecules. These expressions take into account both the molecular 
dielectric anisotropy and the modulation of the local polarisability by density fluctuations. They also 
incorporate all the retardation effects which occur in such liquids. We show that the spectrum splits 
into a g-independent rotational contribution and g-dependent term, which reflects the propagation 
of longitudinal phonons. In the latter, the two light scattering mechanisms enter on an equal 
footing and generate three scattering channels. We study the influence of the two new channels and 
show that they may substantially modify the Brillouin line-shape when the relaxation time of the 
supercooled liquid and the phonon period are of the same order of magnitude. 

PACS numbers: 64.70 Pf Glass transitions - 78.35.+C Brillouin and Raylcigh scattering; other light scattering 
— 61. 25. Em Molecular Liquids 



I. INTRODUCTION 

The theory of inelastic, low frequency, light scattering 
in molecular liquids has gone through different episodes 
over the years. The presence of a narrow inelastic dou- 
blet in the light scattering spectrum due to propagating 
sound waves was predicted by Brillouin in 1922 Its 
first observation was reported in 1930 by Gross, who no- 
ticed that the spectrum consisted of the Brillouin doublet 
superimposed on a central Rayleigh line [|| . An explana- 
tion of this Rayleigh-Brillouin triplet, based on the equa- 
tions of macroscopic hydrodynamics, was suggested by 
Landau and Placzek in 1934 ||^ and was given in com- 
plete form by Mountain in 1966 Q. This hydrodynamic 
description did not take into account the nature of the 
constituents of the liquid. Yet, it can play an important 
role and typical molecular effects have been detected on 
those low frequency spectra (^50 GHz). Let us briefly 
mention three of them. 

One takes place with molecules that are non-spherical 
tops. In that case, both the molecular polarisability ten- 
sor and the tensor of inertia are anisotropic and they can 
reasonably be considered to depend only on the molec- 
ular orientation. The (collective) orientational dynamics 
of the molecules can thus be detected through the cor- 
responding fluctuations of the anisotropic part of the di- 
electric tensor. This corresponds to a central peak with a 
line shape independent both of the scattering vector and 
of the polarisation of the incident and scattered beams. 
Conversely, the tensorial character of the detection mech- 
anism (anisotropic part of the local dielectric tensor) im- 
poses specific relationships between the intensities of the 
spectra when recorded in different geometries. In partic- 
ular, when, as is usually the case, dipole induced dipole 
effects (DID) can be neglected ||], /yy(tj) — ^Ivh{<^)- 
Here the first index corresponds to a vertical (V) polar- 
isation of the incident beam, and the second index to a 



vertical (V) or horizontal (H) polarisation of the scat- 
tered beam, the scattering plane being horizontal. 

The second effect is the possible coupling of the density 
fluctuations with some internal degree of freedom of the 
molecule, having a temperature dependent lifetime, t{T). 
Mountain Q showed that, when ujbt{T)^ 1, where ujb 
is the longitudinal phonon frequency [ p3[ , an additional 
central peak appears in the VV spectrum, which can be 
formally interpreted as the introduction of a retardation 
effect in the bulk viscosity of the liquid. 

A third effect showed up in the study of the VH spec- 
trum of viscous molecular liquids formed of anisotropic 
molecules. A wave-vector dependent feature was de- 
tected, which altered the line-shape of the central peak 
at moderate viscosities |^ and transformed into an un- 
derdamped transverse phonon spectrum for higher vis- 
cosities . It was rapidly recognised that the additional 
spectrum was, indeed, due to a transverse, diffusive Q , or 
propagative Q excitation: any local shear motion inside 
the liquid couples to the mean local molecular orienta- 
tion, and renders the latter anisotropic. This corresponds 
to a similar change of the mean molecular polarisability 
tensor. The detection mechanism of these shear modes 
is thus, as for the first effect, a change in the molecular 
orientations. 

To rationalise the results obtained in VH light scatter- 
ing experiments, some of the present authors (C.D. and 
R.M.P) and co-workers, [l^ here after referred to as 
[I] , proposed to somewhat generalise the usual hydrody- 
namic equations by incorporating through phenomeno- 
logical arguments, in a systematic way: 

a) the coupling of the shear deformation to a local 
mean orientation of the molecules; 

b) a retardation effect in each term corresponding to 
the damping of a variable. 

Concentrating on the case of axially symmetric 
molecules, they characterised each of them by the ori- 
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entation of its axis u, with polar angles 6, 4>, represented 
by P{0, (p, r, t) the local probability density of finding u in 
that direction, and defined a set of orientational density 
variables, Qij{r,t) by: 

Qy(r,i) = J sm0d0d(PP{0,(P,r,t)[u,u, - i<5y] , (1) 

where latin indices i, j, .. represent cartesian components. 
The set of Qij{r, t) forms a symmetrical, traceless second 
rank tensor. 

It was proposed in [I] that the hydrodynamic equations 
pertinent to the case of a such a molecular viscous liquid 
would consist, after linearisation, of the two conservation 
laws: 

p{r,t) + dkJk{v,t) ^0, (2a) 



jk{r,t) = diaki{r,t) , (2b) 

where p(r,t) is the mass density, J{r,t) is the mass cur- 
rent density. Furthermore, they suggested the constitu- 
tive equations for the stress tensor: 

(Tij = {-6P -f T^f, (g) dkVk)Sij +1]^® Tij - fii^ Qij , (3) 

and for the orientation: 

Q,j = -uj^Qrj - r' (g) Qy- + AV ® . (4) 

Here the momentum density J is related to the velocity 
field V via the mean mass density pm- 

Ji{Y,t) = PmVi{T,t) . (5) 

The strain rate (r, t) is a second rank symmetric and 
traceless tensor defined locally by: 

2 

Tij = djVi + diVj - -SijdkVk ■ (6) 

The pressure change 6P{r, t) is related to the instanta- 
neous mass density change Sp(r, t) by: 

SPir,t)^c^Sp{v,t), (7) 

where c is the relaxed sound velocity. The retarded 
couplings are given in terms of ?7f,,rys,/x and T' and 
are, respectively, the bulk and shear viscosities, the 
rotation-translation coupling and the orientational relax- 
ation functions, the symbol (g standing for a convolution 
product in time. Finally, wq is the libration frequency 
of the axial molecules and A' is the rotation-translation 
coupling constant, a quantity that takes into account the 
fact that p and Qij have different dimensions. 

Noting that, for motions that do not involve a local 
density change, the local dielectric tensor can be written, 
neglecting DID effects, as: 

5ey(r,i) =6gy(r,t), (8) 



it was shown in [I] that the preceeding set of equations 
leads to an expression of the intensity that describes the 
complete thermal evolution of the VH spectrum of a su- 
percooled molecular liquid. 

Recently, another of the present authors (A.L.) and co- 
worker made an additional step. Within the frame- 
work of the Molecular Mode Coupling Theory (MMCT) 
|]l2| , they expressed the complete dynamics of a system 
of linear molecules characterised by the position of their 
center of mass and their orientation. They showed that, 
in the q ^ limit, both the correlation functions of 
the density fluctuations and of the orientation fluctua- 
tions contributed to the light scattering mechanism in 
the VV geometry. This result stressed the existence of a 
coupling between these two variables for a longitudinal 
phonon. Nevertheless, the expression they obtained for 
the VV intensity did not make clear two aspects. One 
was the separation of the intensity into a g-independent 
term, representing a pure rotational dynamics, and a hy- 
drodynamic contribution. The second was that the lat- 
ter could be factorized into two parts: the longitudinal 
phonon propagator and a phonon-photon scattering term 
in which both the density and the orientational fluctua- 
tions enter on an equal footing. Conversely, in another 
paper published soon after by two of the present authors 
(T.F. and A.L.) and co-workers this separation into 
two different contributions and the emergence of several 
scattering channels was made apparent. Yet, the tech- 
nique used in ]l3| did not imply any specific light scatter- 
ing mechanism; it was thus not possible to predict from 
that paper which terms would dominate the spectrum for 
a given physical system. 

The present series of papers aims at completing the 
picture of low frequency scattering in viscous molecular 
liquids through the combination of two different aspects. 
On the one hand, we shall present a full length, first prin- 
ciple derivation of Eqs. (||J^) through a Zwanzig-Mori ap- 
proach and obtain from this technique the corresponding 
Onsager conditions which, when met by the relaxation 
functions, guarantee the Brillouin spectra to have a pos- 
itive value whatever the frequency. On the other hand, 
we shall complete the results obtained in [I] by deriving, 
using the same phenomenological equations, expressions 
for the intensities which can be measured in the VV and 
HH scattering geometries; those are the geometries gen- 
erally used to study longitudinal phonons. The results 
will be completely expressed with the help of the differ- 
ent quantities entering Eqs. and a generalisation of 
Eq. (^) which includes density fluctuations. 

This second series of results have immediate implica- 
tions for the experimentalist and they can be derived 
through elementary algebraic techniques. The formal 
proof of the validity of the phenomenological equations 
and of the conditions under which the spectra are positive 
requires the use of more elaborate tools. In order to make 
the present results accessible to as large an audience as 
possible, we found it useful to reverse the logical order 
and to split our presentation into two consecutive papers. 
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The second one (Part II) will give a complete derivation 
of these equations and of their consequences, only briefly 
sketched in , as well as a comparison between the re- 
sults obtained through the phenomenological approach 
and through the more abstract technique of [|l3j. The 
phenomenological part of our work (Part I) is organised 
as follows. 

In Section II, we will derive the expressions for the in- 
tensity obtained in VV and HH experiments. Section III 
will discuss the changes in the line-shape of the longitu- 
dinal phonons that can be expected when the coupling of 
the rotation of the molecules to the longitudinal phonons 
is taken into account in the expression of the dielectric 
fluctuation. Finally, a brief summary of these results, 
a comparison with the expressions previously obtained 
for the VH intensity [I] and additional remarks will be 
presented in Section IV. 



II. THE LONGITUDINAL PHONON 
BRILLOUIN SCATTERING PROBLEM 



In a light scattering experiment, the incident laser may 
be characterised by the amplitude of the electric field Ei , 
and its polarization e^. The spatial and temporal varia- 
tion of the electric field Ej (r , t ) = EiBi exp i (kj • r — o^ii ) is 
given in terms of the wave vector ki and the frequency Ui. 
From the corresponding quantities of the scattered beam, 
et,kt, Wf, one obtains insight into the fluctuations of the 
sample that occur at the scattering vector q = kj — kf at 
the frequency shift uj = LUi — ojf. The thermal fluctuations 
of the dielectric tensor 5eij{r,t) can be decomposed into 
its spatial Fourier components: 



Seij (q, t) ^ (fr de^j (r, t) exp(iq • r) , 



(9) 



a notation that we will use also for other quantities for 
the rest of this paper. This dielectric modulation rep- 
resents a momentary grating by which the laser light is 
scattered. Due to the polarisation of the incident beam 
and of the analyser for the scattered one, the detector col- 
lects only fluctuations corresponding to the projection of 
the dielectric fluctuations onto the two polarisations: 



(5efi(q, i) = effc(5efc;(q, <)ei, 



(10) 



Since the fluctuations readily disappear, there is a cor- 
responding frequency shift, w, leading to a total scattered 
intensity: 



-^fi(q,w) 



dt{Sefi{q,t)Sel{qy)cosiL0t), (11) 



with the notation Se'^ (q) — Scfi (q, t = 0) and similarly 
for other quantities. Since in real space the fluctuations 
are real, one finds (5e^(q)* — (5e^(— q). Furthermore we 
left out well-known factors that can be found, e.g. in the 
book of Berne and Pecora . 



In the case of longitudinal phonons, the expression for 
Scij used in Eq. (||) is incomplete. A contribution pro- 
portional to the mass density fluctuation, Sp, has to be 
added, leading to: 



Seij{r, t) = a5p{r, t)8.i.j + bQij{r, t) , 



(12) 



where Sp and Qij are here the density change and its 
molecular orientation counterpart. Their respective spa- 
tial Fourier transforms couple dielectric fluctuations to 
longitudinal phonons. 

Since the constitutive equations for the stress tensor 
and the equation of motion of the orientation are in the 
form of integro-differential equations it is convenient to 
use the Laplace- Transform (LT), which we use with the 
convention for functions f{t): 



LT[f(t)]{uj) = i / dtf{t) exp(-iwi) . 



(13) 



The scattered intensity 1^ (q, u) can thus be extracted 
as the imaginary part of: 

X(ei, ef , q, lo) ^ LT[(Seti{q, <)Je° (q)*)] (c^) , (14) 

which defines the light-scattering problem in terms of 
correlation functions. 



A. Expressions for the VV geometry 

In the rest of this paper, we shall make use of the usual 
Berne and Pecora axes for light scattering, renamed, in 
agreement with [I], || for the direction of q, _L for the 
direction perpendicular to the scattering plane, and _L' 
for the direction perpendicular to || and _L. We are in- 
terested here in the usual VV scattering geometry where 
Gi and Gf are parallel to _L so that the r.h.s. of Eq. (|l0| ) 
reduces to a6p + bQ±±. The expression of x(ei, ef , q, u) 
is thus completely determined by the knowledge of the 
Laplace transform of the four correlation functions: 

(<5p(q,<)<5pO(q)*),(Mq,i)Q%(q)*),(Q±±(q,i)^p°(q)*> 

and (Q±_L(q, t)(3j__L(q)*), that we shall compute in turn. 

The spatial Fourier transforms of Equations (2a) and 
(Eq) can be grouped into the single equation: 



'5p(q,0 = gVii ||(q,t) . 



(15) 



Inserting the constitutive equation for cr||||, Eq. (^), and 
performing the Laplace transform yields with the help of 
Eqs. d) and see below: 



c^q^Sp{q,uj) - q^^i{uj) wQ|||| (q, w) - g|||| (q) = 

\iodpiq,u;)-Sp°iq)]{i6) 



q- 4 
— Vlbi^) + i;ris{uj)\ - w 

Pm 'J 



Here, we left out the term containing fP{c\) since it will 
drop out once correlation functions with variables of even 
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time parity are built. Due to the coupling we also need 
the Fourier-Laplace transform of Eq. ^ : 

[T'iu)-u] [c.Q,,(q,c.)-Q^/q)] . (17) 

As above, the term containing Q", (q) has been dropped 
since we will correlate Eq. ( ]l7| ) with quantities of even 
time parity only. Using: 
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[c.<5p(q,c.)- V(q)] , (18) 



Tj_^(q,cj) = -igw||(q,t^) 



2i 



[c.Mq,^)-'5p"(q)] , (19) 



one obtains from Eq. (|T^) for the orientational fluctua- 
tions: 



Qiiii(q,^) 



4A' 

'iPm 



r{cu) [<5p(q,c.)-<5pO(q)H 



D(w) J uj 



(20) 



2A' 



^j_j_{q,uj) = r{uj)[Sp{q,u) -Sp°{ci)/uj] 

•jPm 

D[ljj) J LU 

We have introduced, here, the quantity D{uj) which de- 
termines the frequency dependence of the pure orienta- 
tional motions: 



D{ijj) = 0^5 + bjT' [uj) — UJ 



(22a) 



and for later use we also introduce the rotation- 
translation coupling: 



r(w) = uJii{ijj)[D{uj)]-^ 



(22b) 



From Eqs. ( |l6| , pO| , plD one can solve for orientational and 
density fluctuations in terms of their respective initial 
values: 

u5p{q,Lo) = gM(q,c^)-H5pO(q), (23) 



t^OiiiKq-t^) 



D{uj) 



QUl|(q) 



4A'r(w) . 

opm 



(24) 



w(9±j_(q,cj) 



2AV(tj) , 

-—^q^A{q,u), (25) 



with: 



A{q,uj) = Pdq^uo) \c'5p'\q) - Wo^r(cc>)Q|[|| (q)J . (26) 

Here PL{q,uj) abbreviates the longitudinal phonon prop- 
agator: 



PL{q,Lo) = [c^^ - q^c^ - q^ujriL{uj)/ Pji 



(27) 



and, in this equation, the coupling to the parallel compo- 
nent of all the damping mechanisms is expressed in terms 
of the longitudinal viscosity: 



?7L(tj) = 776(cj) + - 



T^s{^)-—D{uj)r(ujf 

LO 



(28) 



The density-density and the density-orientation correla- 
tion functions are then obtained from Eqs. ( |2^j2^ ): 

LT[{Sp{q,t)Sp'{qr)KLu) = 

-[l + q'c'PL{q,ij)] (|p°(q)P), (29a) 



LT[{Sp{q,t)QlAqr)]iq,Lu) = 
_^2ii!Mp^(q,^)(Q0|(q)*Q0^^(q)) , (29b) 



once one has taken into account that both (|p''(q)p) and 
(*9|||| (q)*'3j_i(q)) are of order unity, while: 

(5p°(q)*Q%(q)) = (g%(q)*5p°(q)) = 0{qH') , (30) 

where I is a typical intermolecular distance. The preceed- 
ing relation derives from the fact that one can write, for 
instance: 



(Ql|||(q)V°(q)) = 

N N 
a=l 13=1 



q -2 



(31) 



where Rq. (resp. R^g) are molecular centre of mass posi- 
tions and M^ii is the projection of the unit vector u of the 
molecule /3 on the direction ||. The exponential factor 
in the preceeding equation can be expanded in powers 
of q • (Rq — R^))- One easily convinces oneself that, for 
symmetry reasons, the first two terms of the expansion 
average to zero so that the first nontrivial contribution 
comes from the [q • (R^ — R/3)]^ term, which is 0{q^). 
Conversely, using the same argument for the averages 
appearing in Eqs. (p9|), one finds that the first term 
of the expansion already gives a nonzero value in the 
long- wavelength limit and is proportional to the absolute 
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temperature. The dependence on q can thus be omitted 
in those expressions and will also be, in the rest of the 
paper, for all those equal-time averages. 

For the VV scattering one has to study also the auto- 
correlation of the J L component in the orientation. 

From Eq. (|25|) one obtains: 



Inserting these results into Eqs. (|Tl|,|T^, one finally ob- 
tains for the intensity in a VV scattering experiment: 



LT[(Q^^(q,t)5pO(q)*)]H = 

2A' r(LU) n 9 „ , ^ ,, n,9, 

3pm 



(32) 



Lr[(Qi^(q,t)g%(q)*)](c.) 



_L_LI 



2A' 



(33) 



Here some comments are in order. 

a) One readily verifies that, would one have taken the 
_L_L' component of Eq. ( pj|) instead of its _L_L component, 
one would have obtained: 



iT[(Qii-(q,<)g% 
1 



- 1 - — ^ 



(q)*>]H 



(34) 



Up to a &^ factor, the imaginary part of the r.h.s. of the 
preceeding equation is simply what is referred to as the 
back scattering depolarised spectrum, Ivh{^), which is 
not sensitive to q for long wavelengths. The imaginary 
part of the first term of Eq. (|3^ ) is thus the non-acoustic 
part Ivv{^)i up to the same factor (see Eq. (p^). 
Furthermore, since Qij is a traceless tensor of order two 



{\QW) = -M^'?) 



(35a) 



the first term of the r.h.s of Eq. (33) represents the well- 
known result that the g-independent part of lyv coin- 
cides with 4/3 of the back scattering spectrum Ivh{'-^)- 
b) The four other terms, Eqs. ( p9| , |3^ ) and the second 
term of the r.h.s of Eq. ( |33| ) all contain the phonon prop- 
agator, Pi(g,a;), and are thus g-dependent. Also, due to 
the spherical symmetry of the liquid: 



(Q^x(q,t)<5p"(q)*) - (^p(q,i)Q%(q)* 



(35b) 



f 4fe2 


1 - 




1 3 







A' 



}{uJoqfPL{q,uj) 



2A' 

3pm 



br{uj) 



(36) 



where the preceeding formula makes clear that the first 
term is 4/3 the Ivh{^) back scattering spectrum. 
The preceeding results call for two remarks. 



a) Equation (35c) simply expresses the equipartition 
of energy between the centre of mass motions and the 
libration motions: 



4A' ^1*^11 



Pm 



oc Kb 



(37) 



and this is a necessary condition for the consistency of 
the phenomenological theory summarised in Section I. 
Equation (^^ will appear in a natural way in the micro- 
scopic derivation of the phenomenological equations (see 
Eqs. (27c,28b) of Part II). 

b) Equation ( p6| ) completes previous results obtained 
in [|ll| , separating clearly the role of the phonon propa- 
gator, Pi^{q,uj) from that of the scattering mechanisms 
(a and 6, with appropriate factors). It also separates the 
non-hydrodynamic, rotational contribution (first term of 
Eq. (|3^)) from the hydrodynamic one. In the latter, 
one part (term in b^) is entirely due to scattering by 
the anisotropic part of the polarisability tensor of the 
molecules. The same mechanism was at the origin of the 
g-dependent part oi lyiiiq,'-^) in [I] and in both cases the 
scattered intensity is proportional to a phonon propaga- 
tor, longitudinal or transverse, multiplied by the same 
r{uj)^ factor; the power two stresses that the molecular 
orientation acts twice, once as the source of the fluc- 
tuation and, the second time, as the detection mecha- 
nism. The more complex form obtained here for the q- 
dependent part of /yy((7, w) is the direct consequence of 
the existence of two parallel channels, dp and Qij, in the 
case of longitudinal phonons; this is in contradistinction 
with the single channel case of the transverse phonons 

[I]- 



Comparing Eqs. (29b) and (|32|), the preceeding relation 
implies: 



{QuQl±) 



Hp' 



(35c) 



while, because Qij is a traceless tensor, one also has the 
relationship: 



B. Results for the HH geoemtry 

The same technique can be used to derive the intensity 
which can be obtained in an HH experiment; indeed in 
such a case: 



{QuQW)^-l{\Qn\') 



(35d) 



SeHH{(l,t) = Se±-±-{q,t) sin^ - - (5e|||| (q, t) cos^ 2 ' 
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where is the scattering angle. Expanding in terms of 
density and orientation, one has to consider: 



~a5p{(\, t) cos 9 + bQj^i^i (q, t) sin^ - 



III. ANALYTICAL AND NUMERICAL 
DISCUSSION 

A. Introduction 



l|(q,i)cos -. 



(39) 



To calculate the corresponding auto-correlation function 
as required by Eq. (pi]), one needs again the dynamics 
of density and orientation. Repeating the calculations of 
Sec. II. A, one convinces oneself that Eq. ( [25| ) remains 
valid if _L_L is replaced by _L'_L'. Due to rotational sym- 
metr y, the cross correlators are again identical as in Eq. 
( |35b| ) . Of the six correlation functions that can be built 
in Eq. (^9|), one can easily compute the missing three 
by the methods of the preceeding subsection. From Eq. 
(P4|), one finds: 



Lr[(Q||||(q,i)QO,^,(q)*)](c.) 



ir[(Qiiii(q,i)Qlii|(q)*>]M = 



1 - 



-— ^r(a;)2g2p^(q,^) 

opm 



(40) 



(41) 



In this Section, we discuss some typical features of the 
spectra obtained in a VV scattering experiment. As we 
shall see, the rotation-translation coupling gives rise to 
spectra which differ in shape and, more importantly, in 
some cases in the values for the relaxation times, would 
they be extracted from the widely used model where den- 
sity fluctuations arc considered as the sole light scatter- 
ing mechanism for the g-dependent part of the polarized 
spectrum, i.e. : 



jsinqle / \ 



-{\pX)ImPL{q,uj). (44) 



In order to demonstrate the effects which take place 
when translation-rotation coupling is non-negligible, we 
generate data for the g-dependent part of a VV experi- 
ment by using the full expression, Eq. (|3^). The spectral 
shape is then determined by the frequency dependence of 
the four fundamental memory kernels r]s{uj),r]b(uj), ii(uj) 
and T'{uj). To keep the discussion as simple as possible, 
we shall model the relaxation kernels by single exponen- 
tial Debye processes: 



Lr[(Q||||(q,t)<5pO(q)*>](c.) = 



. , (42) 

Collecting all the terms appearing in Eqs. ( p^ , |39| ) with 
Scfi = SeHH, one arrives at our final expression for the 
HH scattering spectrum: 



, Pm^o 2 n r \ 




,2 1 



1 - 



a cos 6 



bK'r{uj) 

3pm 



(3 4- cos ( 



.(43) 



For a given q, all other scattering geometries give in- 
tensities, which are linear combinations of /yy(g,w) and 
Ihh{<1,^) as well as of Ivniq,^), whose g-dependent 
part detects the transverse phonon, whereas the q- 
independent part detects again lynii^), Eq. (|3^). It 
is thus not possible to detect independently the contri- 
butions of the three terms entering, e.g. Eq. (|3^). Equa- 
tions (|3^) and ( ^ ) are very compact expressions in which 
the g-dependent and the g-independent parts have been 
separated. They allow to study, within some approxima- 
tions for the different quantities entering them, the influ- 
ence of the rotation-translation coupling on the shape of 
a typical longitudinal phonon spectrum. Using for those 
quantities the same toy model as in Jill , we shall proceed 
to this comparison in the next Section. 



M{uj) 



1 



lUJT 



(45) 



where M is any of rig^rji,, fj, or F'. Once the spectra have 
been generated, we treat them as raw data and analyse 
them by a 'density-only' model, Eq. (Q), in which, see 
Eq. (27), ujl = cq and r]i^{ijj) have to be fitted; in the 
spirit of the simplifying assumptions made above, riL{uj) 
is also a simple Debye process, characterised as in Eq. 
(^), by an amplitude, A|^, and a fitted relaxation time, 

TL- 

Let us first make some comments on the general shape 
of the spectra. For a given wave vector q, the second 
part of the r.h.s of Eqs. (^ or ^) is the sum of three 
terms that correspond to three channels of detection of 
the longitudinal phonons: a density-only channel, pro- 
portional to a^, an orientation-only one (oc 6^) and a 
cross channel of density and orientation (oc ab). For a 
simple liquid, the first one is predominant, which depends 
solely on PL(q,u}) and the shape of the corresponding 
spectrum has been studied in detail, in particular for su- 
percooled liquids, in many papers jl^, |l^, 17, |l^, 
The purpose of these papers was basically to extract 
information on the memory function 77^(0;), and on its 
dependence on temperature. In particular, if r]L{oj) is 
characterised by a single Debye relaxation process with 
a relaxation time and an amplitude such that, what- 
ever T^, /mPi(g,u;)/w exhibits a well-defined peak, the 
density spectrum varies in a specific manner with r^: 



7 



for lobTl *C 1, with lob = cq, the spectrum contains a 
quite narrow Brillouin peak centered around, say ujbi, 
with a weak background for other frequencies. The same 
is true for ojbTl ^ 1, but the peak is now centered at 
tOB2 2> iUBi- The only regime where a rather precise 
information on tl is obtained occurs for lubtl ~ 1, i-e. 
approximately 1/5^; lubTl^ 5. 



The influence of the orientation-only channel on the 
shape of the o-dependent part of the VV spectrum was 
discussed in Q , with the help of a similarly simple model 
for the four relaxation kernels; conversely, the cross chan- 
nel ((X ab) was not studied in that paper. We will thus 
first compare the shape of the contributions of the three 
channels. We shall then discuss the line-shape of the total 
g-dependent part of the VV spectrum, focusing mostly on 
the lubt ^ 1 regime. To remain as close as possible to the 
results presented in , we take the same numerical val- 
ues for the various temperature independent parameters, 
c — 0.6, q — 0.02,070 — 1, the units being chosen such 
that A' = 1 and pm = 1. In order to restrict the number 
of parameters, the relaxation time r is the same for all 
four memory kernels indicated in Eq. (|45|), and we use 
the same values of r as in simply studying also and 
in more details the ujbt ~ 1 regime. The amplitudes are 
taken as A^, = 1/V2, A,,^ = and A^ = _y3/16 in 

order to use the same phonon propagator as in . Since 
the Brillouin peak is at much smaller frequencies than the 
libration frequency, loq, the spectra are not sensitive to 
the value of Ap' in the frequency window studied here. 
Hence, we set Ap' = 0, or r'(a;) = 0. Finally, one needs 
a plausible numerical value for the ratio 2A'a/3pm&- Its 
estimate can be obtained from the following considera- 
tions. Impulsive Stimulated Thermal Scattering (ISTS) 
experiments have been performed on supercooled liquids 
formed of anisotropic molecules such as salol and OTP. 
It has been recognised recently that the spectral shape of 
the signal obtained in this pump-probe experiment is sen- 
sitive to the polarisation of the probe beam |l[. In 
these experiments, the origin of the signal is not the ther- 
mal fluctuations of Sp and Qij but the intensity of the 
pump beam; conversely, the signal is detected through 
the dependence of the local dielectric tensor on density 
and orientation as in a VV experiment. The ISTS sig- 
nal then turns out to be linear in dp and Qij with the 
same coefficients as in Eq. (|3^). Making use of differ- 
ent polarisations of the probe beam, the contribution of 
the two variables can be disentangled and a value 
for a typical molecular liquid, metatoluidinc, of the order 
of 0.5 can be derived from the results of ||22l. In order 
to overemphasise the effect of the additional scattering 
channels, we have found it convenient, in the present pa- 
per, to choose a value equal to 4/-\/27 ~ 0.77. All the 
results presented in this Section have been obtained with 
those numerical values. 



B. Analysis of the role of the two additional 
channels 

1. The uobt ^ 1 and the uobt ^ 1 regimes 

• Preliminary remarks 

As already alluded in part A of this Section, the two 
regimes ojbt and w^r 3> 1 are rather uninteresting 
from the practical point of view of determining a value 
of r: the most important piece of information is con- 
tained in the position of the Brillouin peak, its linewidth 
being much more difficult to analyse because its exact 
line-shape cannot be determined. The remaining low in- 
tensity part of the Brillouin line cannot be studied exper- 
imentally because it cannot be disentangled, either from 
the wing of the central peak (first term of the r.h.s of 
Eq. (^)) for high temperatures, i.e. w^r ^ 1, or from 
a fiat background for low temperature, ujbt 3> 1. Let 
us just mention here that, even in these regimes, at very 
low frequency, the r.h.s of the second term of Eq. ( |3^ ) 
can lead to a negative contribution to the intensity, due 
to the existence of the r(u;) and r{ijjY terms, as was al- 
ready mentioned for the latter case in |Tl|| . This is not an 
artefact of the theory, Eqs. but a consequence 

of the introduction of translation-rotation coupling into 
these equations. If one takes it into account, it is no 
longer necessary that the g-dependent part of Eq. ( |3^ ) 
is always positive, the requirement being that only the 
full r.h.s of this equation is positive. We shall differ the 
discussion of this point to Section IV and, mostly, to Part 
II [||. 

In spite of the preceeding remark on the (5-like shape 
of those spectra in the two regimes, their study turns 
out to be instructive, because the way the two addi- 
tional channels contribute to the spectrum will be sim- 
ilar in the ujbt ^ 1 regime we shall study later; there, 
one will directly detect their influence on the line-shape 
while their analytical study is easier in the extreme 
cases we consider here. We shall first develop ana- 
lytic formulae, then demonstrate that these effects do 
show up in the line-shape of the Im[r{u!)PL[q, lo) / lo] and 
Im[r{uj)^ PL{q, uj)/uj\ spectra, if analysed with sufficiently 
large accuracy, and finally show that, nevertheless, the q- 
dependent part of the VV spectrum is little affected by 
these effects in the two cases, (jJ^t <C 1 and cj^r 3> 1, in 
the region of the Brillouin peak. 

• Analytical results 

Table I displays the analytical form of the leading 
terms for the three channels, viz. pure density {a?), pure 
rotation (5^) and cross channel {ab), both for uj ^ lob 
and a; » cjs; for the sake of completeness, the Table also 
displays the approximate Brillouin peak intensity related 
to the first channel. In order to make this Table as easy 
to read as possible, the quantities for the and ab chan- 
nels are the factors by which the results of the first line 
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LOBT <C; 1 

LO <C LOB Brillouin peak lo 3> lob 


ciJsr >• 1 

<C (j-'s Brillouin peak u ^ cij_b 


a" 


A 2 1 „2 A 2 

A r 1 q /\ T 

2 4 A 2 3 4 

f/ c A q CT uj 


A 2 2 » :i 

A r g A 

222/2|a2\2 2a2 6 


RI 


-{corf -{urf 


1 1 


IR 


•2 2 
\ J ^2 \ ! ^2^2 


A' g'A' 


ah RI 




1 1 


ah IR 


2 2 

C 

A' g^A' 


^2 _|_ ^2 ^2 

A' g'A' 



TABLE I: Analytic form for the different scattering channels in the two Umiting cases. The first line, a^, gives the 
analytic form of Im[PL{q,uj)/uj] for uj <C lob, at the peak value and for uj ^ lub, for the two cases ujbt ^ 1 and ljbt S> 1. The 
last four lines give the factors by which the first line has to be multiplied to obtain the analytical form of the corresponding 
term; the factors are given for the 6^ channel (second and third lines) and for the ab channel (fourth and fifth lines). In both 
cases, the two terms, R(eal) I(maginary) and I(maginary)R(eal), are given in the even, respectively odd, lines (see Eq. (M)). 



(a^ channel) have to be multiplied in order to obtain the 
corresponding contribution to the intensity. 

Furthermore, while the intensity of the pure density 
channel is simply equal to Im[PL{q,uj)/u!], in the case 
of the pure orientation channel for instance, its intensity 
/f,2 (uj) is given by: 

-Im[r{Lu)^PLiq,Lo)] = -\ Re[riu;f]Im[PL{q,Lo)] + 

U) UJ I 

Im[r{ujf]Re[PLiq,Lu)]Y (46) 

and a similar expression holds for the cross channel. Since 
/m[i-L((7, w)/w] has typical functional form [1 + X^]~^ 
with X = UJ — ujb for a; > 0, the corresponding real part 
Re[PL{q,uj)/uj] is described by X/\\ + X^]; the contri- 
butions of the two terms of Eq. (E6|) , and similar terms 
in the ah channel, have different line-shapes and have to 
be considered separately. They are denoted as RI and 
IR contributions in the lower part of Table I. Finally, 
we simplified our study by admitting that, in the lim- 
ited frequency range for which this Table is constructed, 
wb/IO^lj^IOwb, the dimensionless quantities ujt and 
ujbt have always the same order of magnitude. Table 
I shows that one has to study three different cases: 

a) ujbt > 1 

Within the approximation just mentioned, the 
and the ab channels have the same multiplying fac- 
tors, and their RI parts give rise to exactly the 
same line-shape as the channel. For uj ^ ujb, 
the (cj/gA)^ factor of the IR term increases the in- 
tensity with respect to the channel; conversely, 
for UJ <^ ujb, the same IR term gives a negative 
contribution: the relative intensities will be lower 
than in the a^ channel. 

b) UJBT < 1 

- The &^ channel has a general (wr)^ factor 



which renders its contribution small on the whole 
frequency domain whatever its sign. The IR part 
gives a negative contribution for uj < ujb which 
changes into a positive one at some value above 
UJB, while the RI part always remains negative. 

- The ab channel gives a line-shape quite dif- 
ferent from the ones discussed above since it is 
dominated by its IR part with a prefactor of order 
unity while the RI part is suppressed by a factor of 
(cjt)^. This results into a typical X/[l + X"^] line 
shape with a negative contribution for uj < ujb, 
a minimum below ujb, a zero in the vicinity of 
UJB, and a positive contribution above ujb, the 
maximum of the latter being slightly above ujb- 

• Numerical results 

Figures 1,2 and 3 illustrate the preceeding results. Fig- 
ure 1 displays the intensity related to b^ channel, /f,2 , for 
T — 1, 10^, lO'^ and 10^. The crossover regime ujbt ~ 1 
corresponds to r ~ 10^. The different curves exhibit the 
shapes discussed above. Note that for clarity the spec- 
trum corresponding to r = 1 has been enhanced by a 
factor of 100. Yet, the spectra for ujbt ^ 1 exhibit an 
additional, strongly negative feature for uj ^ ujb, which 
was not studied in Table I, as it is related to the ujt <C 1 
part of these spectra. 

The different shapes of the lab spectra are shown in 
Fig. 2. These shapes agree with what can be inferred, 
in the ujbt 3> 1 and ujbt <C 1 regimes, from the results 
shown in Table I. 

Finally, Fig. 3 represents the ratio of the total q- 
dependent part of lyv, Eq. (|36|), called Itot, to the 
channel intensity alone, 1^2 . Here the two spectra are 
normalised to the same intensity at the Brillouin peak. 
The factor (w/gA)^, general to the IR part for uj ujb, 
explains the increase of the ratio in this frequency do- 
main, a signature that will remain also in the ujbt ~ 1 
regime. Also, this ratio always decreases, more or less 



FIG. 3: The ratio of KItot to /„2 for the same values of r as in 
Fig.l. For each spectrum, K is adjusted in a such a way that 
the intensities of the two spectra are identical at the Brillouin 
peak. The horizontal line, drawn for a ratio equal to 1, allows 
to locate these Brillouin peaks for the four cases. 

2. Numerical analysis of the u>bt ~ 1 regime 

In order to look for an experimentally detectable ef- 
fect, we have computed the total g-dependent part of the 
intensity, Itot, for the 10 values of r reported in the first 
column of Table II, all close to the ujbt ~ 1 condition. 
Those 10 spectra are shown in Fig. 4 as dashed lines and 
we have checked that they notably differ from the scaled 
spectra that can be obtained for the same values of t con- 
sidering a line-shape as given by Eq. (^4[). To explore 
more closely the effect of the additional channels, both 
on the line shape and on a possible misinterpretation of 
these spectra in terms of 'density-only' spectra, we have 
tried to fit them with a /m[Pi((7, lo)/lu] expression, where 
the longitudinal viscosity was expressed as: 

q UJr-lL(UJj/Prn = A^—— 1- lUJ-fL ■ (47) 

1 + lUJTL 

A value of 7l, consistent with a previous Brillouin scat- 
tering study of m-toluidine 0, — 10~^, was chosen, 
and fits for the spectra corresponding to the ten different 
values of r were performed, treating Al,lul — cq and tl 
as free parameters. 

The numerical values extracted from the best fits, as 
well as the value of the peak frequency of Itot, ojb, are 
reported in Table II. The fits are displayed in Fig. 4 
and they show that the closer ujbt is to unity, the worse 
is the agreement with the density-only model. Note that 
the value of 2A'a/Spmb has been overestimated to empha- 
size the effect. The discrepancy between the fit function 
and the computed spectra has always the characteristic 
features discussed before: the fit function has a lower in- 
tensity than the computed one above the Brillouin peak 
and a higher one below. 



FIG. 4: Itot (dashed line) for 10 values of r in the vicinity of ojbt = 1, and its fit (full line) with a 1^2 model (see text). The 
values of r and lubt are indicated in each case. 



The trends reported in Table II are worth commenting: 
the ratio r^/r increases with longer relaxation times r, 
i.e. lower temperatures: the apparent relaxation time tl 
is longer than the original one. In Brillouin experiments 
on supercooled liquids, this effect will take place for re- 
laxation times of the order of some nanoseconds, while 
it cannot be identified for much longer relaxation times. 
The efi^ect will then be apparent in the region where the 
curve logr versus l/T has its maximum curvature for 
fragile glass forming liquids. This curvature is the origin 
of the nonzero value of the Vogel-Fulcher temperature: 
an apparent increase of the relaxation times in this time 
window will decrease this curvature and result in an ar- 
tificial decrease of the Vogel-Fulcher temperature. 

Table II also exhibits an effect that has not been re- 
ported in experimental studies of Brillouin spectra of su- 
percooled liquids: a decrease, instead of an increase, of 

H, |l9). This contrasts 



Ai upon cooling |15, 16 



with the familiar trend oi ujl which increases with r as 



is shown in Table II. The unconventional dependence on 
temperature of may be an artefact of the oversimpli- 
fied model for the memory kernels. 



IV. SUMMARY AND FINAL REMARKS 

The present paper aimed at giving complete and trans- 
parent expressions for polarised Brillouin light scattering 
experiments within the framework of the phenomenolog- 
ical equations recalled in the Introduction. These equa- 
tions have been derived through the use of heuristic ar- 
guments in [I] and a demonstration of their validity has 
been briefly sketched in . Their detailed proof will be 
given in the second paper of this series (Part II), where 
the same set of relevant variables, i.e mass density p, ori- 
entation Qij and their respective currents are considered. 

The expressions for the intensity in the VV and HH 
scattering geometry are given by Eqs. (p6|) and (Esl). In 
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r " 


TL ' 






LOB " 


15 


13.4 


0.0210 


0.0125 


0.0121 


20 


19.8 


0.0200 


0.0128 


0.0123 


25 


27.6 


0.0190 


0.0131 


0.0125 


30 


37.0 


0.0180 


0.0133 


0.0132 


50 


67.7 


0.0170 


0.0141 


0.0186 


80 


109 


0.0160 


0.0147 


0.0204 


100 


136 


0.0156 


0.0149 


0.0207 


125 


175 


0.0154 


0.0148 


0.0208 


150 


209 


0.0152 


0.0149 


0.0210 


200 


270 


0.0150 


0.0150 


0.0210 



"initial value 
"fitted value 
'^deduced from Itoti^j) 

TABLE II: Fit parameters of the total spectrum for 10 
different values of r. The first column gives the values of 
T for which the total spectra, hot, are computed. The three 
next columns give the fit values for tl, (see Eq. (^)) and 
ujl, the relaxed phonon frequency, while the last one gives the 
frequency, ojs, of the peak of Itot- Note that, for these values 
of r, Itot exhibits a broad enough maximum (see Fig. 4) for 
ujb to be defined with an uncertainty of order 5 x 10""'. 



both cases, the spectrum naturally sphts into a sum of 
two terms: one describes the pure orientational dynam- 
ics of the molecules, decoupled from density fluctuations; 
the second term involves the propagation of longitudi- 
nal phonons. The corresponding longitudinal viscosity 
comprises contributions from the relaxation of the trans- 
lational and orientational motions. The coupling of these 
two types of motion is characterised by the translation- 
rotation coupling constant. A', and by the frequency- 
dependent rotation-translation function, /i(w), while the 
orientational dynamics is characterised by D(uj) which 
depends very weakly on lo in the region of interest for 
Brillouin scattering studies. The role of A' and r{uj) 
in the spectra is twofold: first, their play a role in the 
phonon propagator and, second, they enter in the detec- 
tion mechanism via the factor [a + 2A'br{Lu)/3prn\- 

The role of the molecular polarisability anisotropy in 
detecting both the uncoupled orientational dynamics and 
the transverse, diffusive or propagative modes with wave 
vector q already appeared in [I]. For the sake of com- 
pleteness, we reproduce here the results obtained in that 
paper under a from that allows for an easy comparison 
with Eq. (M): 



verse phonon propagator: 

PT{q,i^) = [w^ - g^p„^w?7T(t^)] ^ 
characterised by the transverse viscosity: 

A' 

77T(t^) = rjs{uj) D{uj)r{ujY . 



(49) 



(50) 



The absence, in the second square bracket of Eq. ( |4§| ) 
of the factor 2/3 is in agreement with a result of jl^, as 
will be discussed in Part II. 

We studied, in Section III, the change in the spectral 
line-shape brought by taking into account the additional 
light scattering channels. We concentrated on the second 
term of the r.h.s. of Eq. (|3^), mostly discussing the type 
of distortion produced by these new terms in the spectral 
shape in the vicinity of the Brillouin peak. Let us add 
some remarks. 

Firstly, in Section III, the numerical calculation of the 
total intensity was performed using a positive ratio hja 
and we assumed through Eq. ( [44)) that /i(t) could be 
characterised by some amplitude, A^, and a smooth, pos- 
itive, decreasing function of time that approaches zero for 
long times. Contrary to the coefficient a in Eq. (^2|), the 
numbers h and are not always positive, but we want to 
point out that their product 6A^ is commonly a positive 
quantity, for prolate as well as for oblate axial molecules. 
Indeed, on the one hand, the polarisability anisotropy, 
6, is usually positive for prolate molecules giving rise to 
glass-forming liquids and negative for the oblate ones. On 
the other hand, the sign of [lii) is also shape dependent as 
can be inferred from the following argument. Subjected 
to a steady shear flow, , a molecular liquid builds up a 
nonzero stationary mean orientation Qij which is easily 
deduced from Eq. (0): 



1 



(51) 



If, for instance, this flow is in the x direction with a pos- 
itive gradient in the z direction, i.e. Txz > 0, one easily 
convinces oneself that a long axis of the molecule will 
be, on average, parallel to x -|- z, and a short axis paral- 
lel to X — z. Hence, cigar-shaped molecules will exhibit 
Qxz > 0, whereas disk-shaped ones lead to Qxz < 0. 
In other words, the shear flow exerts a torque on the 
molecules so that the sign of: 



^{t)dt = Imfi{LU ~ 0) , 



(52) 



/yff((7, w) 



+ ^(c.„g)^cos^-PT(g, 







'f 


\ ^0 







A 



br{Lo) 

Pm 



, (48) 



where the g-dependent part is now mediated by the trans- 



which is also the sign of A^ for a smoothly decreasing 
ftmction Imp,{ijj), is the same as the sign of h. 

Secondly, the total VV spectrum is the sum of two 
terms in Eq. (|3^): the first is proportional to the 
backscattering VH spectrum and gives always a positive 
contribution to the intensity for all frequencies. Con- 
versely, at very low frequencies, the second, g-dependent, 
term may give rise to a negative contribution. This is 
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in contrast to the density-only model and may serve as 
a simple test to determine whether or not translation- 
orientation coupling is significant. Since the measured 
spectrum should be positive whatever the frequency, one 
would like to know if the proposed phenomenological 
equations ensure this property and, if necessary, what ad- 
ditional requirements have to be imposed on the memory 
kernels to guarantee this positiveness. The phenomeno- 
logical equations allow to derive the appropriate condi- 
tions. However, these conditions will appear more nat- 
urally and in a transparent way in the microscopic ap- 
proach. Therefore their derivation will be postponed to 
Part II and we merely state here the result: the measur- 
able spectra in the three geometries VV, HH, VH, Eqs. 
(^) (^sj) and (^), are positive for all frequencies pro- 
vided that: 

a) the imaginary part of T' {u) , rib{LLj) and ris{u!) are 
positive for all frequencies; 

b) the imaginary part of the translation-rotation cou- 
pling fulfills: 



[Imr'{uj)][Ifmjs{uj)] - A'[/m/i(w)]2 > 0. 



(53) 



This conditions appear as a generalisation of the On- 
sager relations for the dynamics of a coupled system and 
are, as expected, independent of the a/b ratio. This cor- 
roborates the argument that the positiveness of the spec- 
tra should not depend on the relative strength of the two 
scattering mechanisms. In view of the analytical form of 
the memory kernels as given in Eq. (Q), and with the 
numerical values used in Section III for A^j_ , A^^ , A^ and 
A', Eq. (H) is fulfilled for aU frequencies if a|, > 1/2. 
Nevertheless, in the frequency range considered, the role 
of Ar' is important only in the g-independent part of Eq. 
(p6|) (rotational part of the spectrum): 



1 



-Im 



1 - 



D{u;) 



= —Im 



ojT'iuj) 



(54) 



We have checked that the neglect of Ar' in D{uj), both 
in the denominator of the r.h.s. of Eq. (|5^ ) and in r(tj), 
has virtually no influence: it does not change the inten- 
sity and the shape, either of the pure rotational spectrum, 
or of the g-dependent spectrum discussed in Section III. 
Conversely, this neglect allows for the analytic discus- 
sion of Itot performed in that section. This justifies, a 
posteriori, the simplification made there. 

Finally, we already pointed out in the Introduction 
that effects related to energy conservation, in particular 
to thermal diffusion, are neglected in the present paper. 
Some consequences of the heat diffusion process are well 
known since the pioneer work of Landau and Placzek . 
It gives rise to contributions in the propagator of longi- 
tudinal phonons for very low frequencies. Other aspects 
of the role of temperature fluctuations were taken into 
account in There, it was shown that, if one deals 

explicitly with the coupling of temperature to dielectric 
fluctuations, one arrives at a rather complex form of the 
coupling of detection channels and propagators of these 
excitations. In order to describe this aspect within a 
phenomenological approach, the constitutive equations 
of the present paper have to be generalised and supple- 
mented by an equation of motion for the energy conser- 
vation, while the detection mechanism, Eq. (|T^), will 
remain unchanged. 
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